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(57) Abstract: The present invention relates 
to a method of estimating the properties of a 
polymer product by using converted process 
variables (cPV's) which means process variable 
that final or intermediate product experienced 
earlier in average in the reactors by stage with 
respect to the final product, which incorporates 
residence time distribution to process variables, 
wherein said method comprises the following 
steps of: computing cPV's by incorporating 
the residence time distribution to said process 
variables by means of theoretically determining 
the amounts of content and discharge of the 
product in the process, and then solving the 
. balance equations of a hypothetical substance by 
taking the respective process variables as those 
of the hypothetical substance; and estimating the 
properties of the product by inputting to various 
property estimation models cPV's in the reactors 
by stage with respect to the final product after 
incorporating the residence time distribution to 
said process variables. 
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METHOD OF ESTIMATING THE PROPERTIES OF A POLYMER PRODUCT 

TECHNICAL FIELD 

The present invention relates to a method of estimating the properties of a polymer 
5 product, or more particularly to a method of estimating the properties of a product by 
using converted process variables (cPV's) which means process variable that final or 
intermediate product experienced earlier in average in the reactors by stage with respect 
to the final product, which incorporates residence time distribution to process variables, 
wherein said method comprises the following steps of: incorporating the residence 

10 time distribution to said process variables by means of theoretically detennining the 
amounts of content and discharge of the product in the process, and then solving the 
balance equations of a hypothetical substance by taking the respective process variables 
as those of the hypothetical substance; and estimating the properties of the product by 
inputting to various property estimation models cPV's in the reactors by stage with 

15 respect to the final product after incorporating the residence time distribution to said 
process variables. 

BACKGROUND ART 

Generally, the adjustment of physical properties of polyolefin in a commercial plant 
20 is largely dependent on experiences of a skilled operator. However, because of the 
differences in modes of operations of respective operators, the operations are 
sometimes inconsistent. The modes of operations which are different from operator to 
operator affect the physical properties of a product (i.e., uniformity) and in the end, 
even if the products are of the same grade, there still may be differences depending on 
25 the lots. Obviously, this is a potential cause of claims or complaints from the end 
users. 

In order to overcome these problems, it is useful to apply the technology, which 
allows estimation of properties of a product, at the point of its production. 

Conventionally, in estimating the properties of a polymer product, various 
30 empirical and statistical models have been mainly used, such as empirical correlation, 
neural network, or partial least square (PLS) models. 
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While the correlation models of prior art can be relatively easily applied to the 
steady states wherein the process conditions are stably maintained, the conventional 
technology is problematic in that it frequently shows large differences from the actual 
analysis values in the unsteady states or in a grade change, where the process conditions 
5 are in the state of flux. 

These differences are due to the existence of a process delay. The process delays 
are expressed in different forms, depending on the intrinsic residence time distribution 
of various devices. For example, there is a simple case of a push back to a certain 
period of time (e.g., in pipes), or there may be a push back of extended effects in the 
10 form an exponential function, such as in continuous stirred tank reactors (CSTR). For 
other reaction devices, there may be other peculiar types of residence time distribution. 

In particular, in cases of multi-stage continuous polymer plants, up until to the 
point that the product would be affected by the extended effects caused by a change in 
operational conditions of the previous stage, there would be. a long time delay, 
15 depending on the intrinsic residence time distribution characteristics of the process. 
Consequently, because of a rather simplistic way of estimating the properties of a 
polymer product from the process variables of the current states, it would be of course 
accompanied by a significant degree of error. 

Meanwhile, a typical method of estimating the properties of polymers by 
20 incorporating the residence time distribution of the process is by means of a physical 
model. It is a method of establishing balance equations for a substance by way of 
reaction system and then solving said equations. For this method based on said 
physical model, it is necessary to obtain reaction rate constants and various types of 
physiochemical constants. However, there lies the problem since it is not so easy to 
25 obtain these types of constants. To solve this problem of physical model, the method 
widely used in the industry is a method of computing cumulative properties (hereinafter 
cumulative properties model) by calculating the instantaneous properties of a product, 
followed by the application of residence time distribution and the mixing rule of 
polymer properties. 

30 The method according to said cumulative properties model has many advantages, 
such as easy application to the actual process in some cases, and potential utilization of 
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empirical correlation, neural network, PLS models, etc. in estimation of instantaneous 
properties. 

However, with respect to the method of cumulative properties model, there is a 
problem of requiring "model training," or the data of the steady states for the purposes 
5 of carrying out the process of empirical correlation, neural network, or PLS for 
estimating the properties of a product. In particular, in case of a system with recycle 
streams which are not maintained constantly by time, it is almost impossible to apply 
the cumulative properties model. 

10 DISCLOSURE OF THE INVENTION 

In order to solve the problems of prior art as above, the present invention seeks to 
provide a method of applying cPV's in operation in the reactors by stage, such as 
temperature and pressure, to the models of estimating the properties of a product by 
processing the process variables by incorporating the residence time distribution to the 

15 process variables themselves by way of further simplification of the cumulative 
properties model which finds cumulative properties by means of residence time 
distribution and the mixing rule of polymer properties by computing the instantaneous 
properties from the current values of process variables. 

The present invention relates to a method of estimating the properties of a product 

20 by using cPV's in the reactors by stage with respect to the final product, which 
incorporates residence time distribution to process variables, wherein said method 
comprises the following steps of: computing cPV's in the reactors by stage with 
respect to the final product by incorporating the residence time distribution to said 
process variables by means of theoretically determining the amounts of content and 

25 discharge of the product in the process, and then solving the balance equations of a 
hypothetical substance by taking the respective process variables as those of the 
hypothetical substance; and estimating the properties of the product by inputting to 
various property estimation models cPV's in the reactors by stage with respect to the 
final product after incorporating the residence time distribution to said process variables. 

30 Moreover, in relation to the method of monitoring the state of process incorporating 
the residence time distribution to the process variables, the present invention relates to a 
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method of monitoring the state of process, wherein said method comprises theoretically 
determining the amounts of content and discharge of the product in the process, and 
then incorporating the residence time distribution to said process variables by way of 
the method of solving the balance equations of a hypothetical substance by taking the 
N 5 process variables as those of the hypothetical substance; and monitoring the state of 
process by inputting to various process-state-monitoring models cPV's in the reactors 
by stage with respect to the final product after incorporating the residence time 
distribution to said process variables. 

10 BRIEF DESCRIPTION OF DRAWINGS 

Fig. 1 is a bloc diagram, which illustrates the method of estimating the properties of 
a polymer product according to the present invention. 

Fig. 2 is a graph, which conceptually explains the need for incorporating residence 
time distribution to various process variables according to the present invention, with 
15 an example of a continuous stirred tank reactor (CSTR). 

Fig. 3 is a graph, which shows the raw data of reactor temperatures and pressures as 
a function of time. The reactor temperature is one of the important elements in 
operating the reaction process of high-density polyethylene (HDPE). 

Fig. 4 is a graph, which shows the raw data of mass flow rates of ethylene and 
20 hexane into the inner chamber of a reactor as a function of time. The rate is one of the 
important elements in operating the reaction process of high-density polyethylene 
(HDPE). 

Fig. 5 is a graph, which shows the raw data of mole ratios of hydrogen and ethylene 
(H2/C2) and cPV's as a function of time. The mole ratio is one of the important 
25 elements in operating the reaction process of high-density polyethylene (HDPE). 

Fig. 6 is a graph, which shows the estimated values of melt index of an 
intermediate product as a function of time from the outlet of a first-stage reactor in the 
reaction process of HDPE according to the present invention. 

Fig. 7 is a graph, which shows the estimated values of melt index of the final . 
30 product as a function of time from the outlet of a second-stage reactor in the reaction 
process of HDPE according to the present invention 
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BEST MODE FOR CARRYING OUT THE INVENTION 

The present invention is described in further detail with an example as below: 
Nevertheless, it is described here for illustrative purposes only and does not limit the 
5 present invention in any manner. 

Fig. 1 is a bloc diagram, which shows the method of estimating the properties of a 
polymer product according to the present invention. Fig. 1 shows the process of 
computing cPV's in the reactors by stage with respect to the final product discharged at 
a certain time period with a given operational condition of each stage. This 
10 computation was made with due consideration of residence time distribution. 

The process of computation shown in Fig. 1 is described as below, which computes, 
with due consideration of residence time distribution, cPV's in the reactors by stage 
with respect to the final product discharged at a certain time period. 

First, the raw data such as reactor temperature, pressure, and concentration were 
15 read through a controller, etc. From the raw data read in as such, the balance 
equations (Mathematical Formulas 1-2 as below) of a substance with respect to the 
polymers were constructed. From these balance equations of the substance, M,y could 
be obtained, which is the mass (kg) of the polymers, polymerized in the 1 th reactor and 
passed over thereto, out of the polymers in the j* reactor. Next, W/jj was obtained, 
20 which is the mass flow rate (kg/hr) of the polymers, polymerized in the i & reactor, out 
of the polymers fluxed in into the j* reactor. Then, W 0jJ was obtained; which is the 
mass flow rate (kg/hr) of the polymers, polymerized in the i* reactor, out of the 
polymers discharged from the j* reactor. The mathematical formulas used herein are 
as follows: 

25 

Mathematical Formula 1 
dMi/dt^Wfjj-Wojjdtj) 



30 



Mathematical Formula 2 
dMa = Wi-W oM (i=j) 
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Mathematical Formula 3 

W/Jj = Wow 

Mathematical Formula 4 

Wo,ij = wfijW 0 j 

Mathematical Formula 5 
wfij = Mi/M/i^) 

Mathematical Formula 6 



Mathematical Formula 7 

In the above mathematical formulas, the following conventions are used: 
Mi : Mass (kg) of polymers in the i** 1 reactor; 
Wf ti : Mass flow rate (kg/hr) of polymers of inflow to the i* reactor; 
20 W 0 j : Mass flow rate (kg/hr) of polymers discharged from the i m reactor; 

w/jj : Ratio (kg/kg) of the amounts of polymers, polymerized in the i* reactor, 

out of the polymers in the j* reactor ; 
Mij : Mass (kg) of polymers, polymerized in the i* reactor and passed over 
thereto, out of the polymers in the j* reactor ; 
25 t :Time(hr); 

Wi : Polymerization rate (kg/hr) within the i m reactor; 

W M : Mass flow rate (kg/hr) of polymers, polymerized in the i* reactor, out of the 

polymers fluxed in to the j* reactor; and 
W 0 ,ij : Mass flow rate (kg/hr) of polymers, polymerized in the i m reactor, out of the 
30 polymers discharged from the j* reactor. 
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The above formulas are rearranged as follows: 

Mathematical Formula 8 
5 dMy/dt = M u y6j. 1 -M i /dj(i#) 

Mathematical Formula 9 
dM i j/dt = W i -M i /6 i (i=j) 

10 In the above mathematical formulas, the following conventions are used: 
0i : Average residence time (hr) in the i m reactor; 
My : Mass (kg) of polymers, polymerized in the i* reactor and passed over 

thereto, out of the polymers in the j ^ reactor ; and 
Wi : Polymerization rate (kg/hr) within the i* reactor. 

15 

Various cPV's were computed by applying the numerical values as obtained above 
to the balance equations of the hypothetical substance (taking the process variables as 
those of the hypothetical substance) of various process variables expressed as 
differential equations. To obtain the properties of a product, various cPV's as 
20 computed above were inputted to the empirical correlation, neural network, partial least 
square models, etc. The process as such was repeated for each unit time. 

Moreover, if necessary, the lot average value, lot maximum value, and lot minimum 
value could be obtained by using the numerical values as computed by using cPV's. 
The present invention is now described in further detail with examples with reactor 
25 temperature, which is a typical operational condition therein. 

By taking the reactor temperature as that of the hypothetical substance, the balance 
equations of the hypothetical substance as follows could be established: 



30 



Mathematical Formula 10 
dMijTi/dt^ WmTw - WojjTtj (i ±j) 
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Mathematical Formula 11 
dMuTi/dt^ WiTi - W oM Tn (i = j) 

In the above mathematical formulas, the following conventions are used: 

5 My : Mass (kg) of polymers, polymerized in the i & reactor and passed over 
thereto, out of the polymers in the reactor ; 
ft : Temperature (°C) in the X th reactor; 

Tg : Of the polymers in the j * reactor (or discharged), the average temperature 
(*C) in the i m reactor, experienced by the portion produced in the 
10 reactor; 

t : Time (hr); 
( W( : Polymerization rate (kg/hr) within the I th reactor; 

Wfjj : Mass flow rate (kg/hr) of polymers, polymerized in the i m reactor, out of 
the polymers fluxed in to the j* reactor; and 
15 W 0tij : Mass flow rate (kg/hr) of polymers, polymerized in the i* reactor, out of 

the polymers discharged from the j* reactor. 

Of the polymers in the j* reactor (or discharged), if T cum>ij were taken to be the 
temperatures experienced in the i* reactor by all of the polymers produced from the 
20 first reactor to the i m reactor, the formulas must be changed as follows: 

Mathematical Formula 12 

*t m Am idt = t w f^, H -Ev^ o' * J) 

k** k=l k=l 



25 Mathematical Foimula 13 

M..T .Jdt = dM.T „ Idt = W_ T J-OTT _W ' 

l eum t U 



«CE" A- /* = dM t T mJl Idt =W fJ T amJ _ u . 1 +W,T, -W oi T a 

*=1 



In the above mathematical formulas, the followings conventions are used: 

Tam.fi : Of the polymers in the jth reactor (or discharged), the temperatures (°C) 
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experienced in the i & reactor by all of the polymers produced from the 
first reactor to the i* reactor; 
Mij. : Mass (kg) of polymers, polymerized in the I th reactor and passed over 
thereto, out of the polymers in the j *** reactor ; 
5 W/j : Mass flow rate (kg/hr) of polymers of inflow to the i* reactor ; 

Mi : Mass (kg) of polymers in the i A reactor ; 

Wf p ij : Mass flow rate (kg/hr) of polymers, polymerized in the i A reactor, out of 

the polymers fluxed in to the f 1 reactor; 
W 0 jj : Mass flow rate (kg/hr) of polymers, polymerized in the i* reactor, out of 
10 the polymers discharged from the reactor; 

W 0 ,i : Mass flow rate (kg/hr) of polymers discharged from the i* reactor; 
W( : Polymerization rate (kg/hr) within the i m reactor; and 
Ti : Temperature (°C) in the i* reactor. 

15 Of course, in case of reaction temperature, among two of the mathematical 

formulas of Mathematical Formulas 10-11 and 12-13, the forms of Mathematical 
Formulas 10-11 should be used to confer some physical meaning therein. Depending 
on the characteristics of process variables, there should be appropriate selections for use 
among Mathematical Formulas 10-11 or 12-13. 

20 After going through the procedures as above, the instantaneous operational 
conditions and cPV's at each stage of the production point could be computed. 

Fig. 2 is a graph (case of having residence time distribution of CSTR), which 
conceptually explains the method incorporating the residence time distribution to the 
respective process variables according to the present invention. In other words, it 

25 conceptually shows the method of selecting cPV's at some instant while working with a 
single reactor with the reactor flow corresponding to a CSTR. The operational 
conditions of a product which is discharged at some instant (/) do not correspond to the 
operational variables at that instant (t). Rather, they correspond to the values of sum 
of the values of operational variables respectively multiplied by a certain weight factor 

30 from that instant to some previous point. The sum of the weight factors is 1, and 
accordingly the area indicated as "B" is 1. In case of a single CSTR, the formula of i- 
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- exp[-(t-At)/d] is given. At is a time interval from t to some previous point in time. B 
is the average residence time in the reactor. Theoretically, At can reach infinity, but it 
is generally considered up to the time interval of 3-4 times of 0. 

The operational conditions of a product discharged at an instant t as above do not 
5 correspond to the values of operational variables at that time f. Rather, they 
correspond to the values of sum of the values of operational variables respectively 
multiplied by a certain weight factor from that instant to some previous point. The 
problem caused by this inconsistency becomes more serious in case of. having long 
residence time within a reactor, or in case of using a multi-stage reactor. To solve this 
10 problem, as shown in Fig. 1, it is necessary to employ a method of setting cPV's in the 
reactors by stage with respect- to the final product at some instant in the multi-stage 
reactor. 

The method of estimating the properties of a polymer product as above can also be 
used as a method of monitoring the state of process of the polymer product. In other 

15 words, in relation to the method of monitoring the state of process incorporating the 
residence time distribution to the process variables, the method of monitoring herein 
involves theoretically determining the amounts of content and discharge of the product 
in the process, and then incorporating the residence time distribution to said process 
variables by way of the method of solving the balance equations of a hypothetical 

20 substance by taking the process variables as those of the hypothetical substance; and 
monitoring the state of process by inputting to various process-state- monitoring models 
cPV's in the reactors by stage with respect to the final product incorporating the 
residence time distribution to said process variables. 

Moreover, by using the method of estimating the properties of a polymer product as 

25 above, it can compute the average operational conditions appropriate for optimized 
operation per lot or batch of the polymer products, and it can also be used in quality 
control of polymer products. 

Moreover, by using the method of estimating the properties of a polymer product as 
'above for the purposes of monitoring the state of processes therein, it can compute the 

30 average operational conditions appropriate for optimized operation per lot or batch of 
the polymer products, and it can also be used in quality control of polymer products. 
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Below, the present invention as above is described in detail with an example. 
However, it is for illustrative purposes only and should not be deemed to limit the 
present invention in any manner. 

Figs. 3-7 are graphs according to one embodiment of the present invention. "The 
5 process therein was a high-density polyethylene reaction process in which two reactors 
were connected serially. 

Figs. 3-5 illustrate the changes in raw data (as a function of time) of various 
important factors during the operation of a high-density polyethylene reaction process. 
Fig. 3 shows the changes in temperature and pressure within the reactor. Fig. 4 shows 

10 the changes in mass flow rates of ethylene and hexane as a function of time. There, 
ethylene was the main raw material flowed in to the reactor with a diluent, and hexane 
was used to control the residence time. Moreover, Fig. 5 respectively illustrates the 
raw data of mole ratios of hydrogen and ethylene (H2/C2) in gas phase within the reactor. 
There, hydrogen was used as an agent for controlling the molecular weights of 

15 polymers. The mole ratios of hydrogen and ethylene in gas phase were analyzed using 
an on-line gas chromatograph. As for the other operational conditions of the reactor, 
the capacities of first-stage and second-stage reactors were respectively 71m 3 , and the 
production rate of polymers was approximately 20 ton/hr. The reactor residence time 
for the first-stage reactor was approximately 2 hours, and approximately 1 hour for the 

20 second-stage reactor. In particular, Fig. 5 also shows cPV's as a function of time, in 
addition to the raw data of mole ratios of hydrogen and ethylene (H 2 /C 2 ). As for the 
computation method of cPV's, Mathematical Formulas 10-11 as presented above were 
used herein. As shown in the figures, in contrast to those of the raw data, the graphs 
of cPV's incorporating the residence time distribution of the reactors at each stage 

25 showed a form of a very smooth curve. By inputting cPV's, not the raw data, to the 
empirical correlation, neural network, PLS models, etc., the properties of intermediate 
and final products could be thus obtained. 

Fig. 6 is a graph, which shows the estimated values of melt index (as a function of 
time) of the intermediate products at the outlet of the first-stage reactor during the 

30 HDPE reaction process according to the present invention. 

Fig. 7 is a graph, which shows the estimated values of melt index (as a function of 
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time) of the final products at the outlet of the second-stage reactor during the HDPE 
reaction process according to the present invention. For reference, it also shows the 
estimated values of properties (i.e., instantaneous MT) of the case in which the raw data 
were directly used instead of cPV's. As shown in the figure, the estimated values of 
5 properties (i.e., cumulative MI) based on cPV's have a much better fit against the actual 
analytical results obtained from the laboratory. 

For reference, in the embodiment, the method of empirical correlation was used as 
a properties estimation model. The method of empirical correlation in the 
embodiment is one that expresses melt index of the polymers in the form of a function 
10 of temperature, pressure, and composition ratios of the substituents in the process. 

According to the method of estimating the properties of a product by using cPV's in 
the reactors at each stage with respect to the final product, based on the present 
invention, it can easily estimate the properties of a polymer product from a continuous- 
type polymerization reactor and also monitor the state of each process therein. 
15 Moreover, according to the present invention, even with a batch or semi-batch type 

polymerization reactor, or in the reaction system of uneven treatment of polymers, the 
method can easily estimate the properties of the final product while allowing 
monitoring of the state of each process therein. 



20 
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CLAIMS 

What is claimed is: 

1. A method of estimating the properties of a polymer product, which comprises the 
5 steps of: computing converted process variables (cPV's) which means process 

variable that final or intermediate product experienced earlier in average in the 
reactors by stage with respect to the final product by incorporating the residence 
time distribution to process variables by means of theoretically determining the 
amounts of content and discharge of the product in the process,- and then solving 
10 the balance equations of a hypothetical substance by taking the respective process 

variables as those of the hypothetical substance; and' estimating the properties of 
the product by inputting to various property estimation models cPV's in the 
reactors by stage with respect to the final product after incorporating the residence 
time distribution to said process variables. 

15 

2. The method of estimating the properties of a polymer product according to Claim 
1, wherein said process is a process which uses a single polymer reactor or a 
multi-stage polymer reactor. 

20 3. A method of monitoring the state of process, incorporating the' residence time 
distribution to process variables, wherein said method comprises theoretically 
determining the amounts of content and discharge of the product in the process, 
arid then incorporating the residence time distribution to said process variables by 
way of the method of solving the balance equations of the hypothetical substance 

25 by taking the respective process variables as those of the hypothetical substance; 

and monitoring the state of process therein by inputting to various process-state- 
monitoring models cPV's in the reactors by stage with respect to the final product 
after incorporating the residence time distribution to said process variables. 

30 4. The method of monitoring the state of process according to Claim 3, wherein said 
process is a process which uses a single polymer reactor or a multi-stage polymer 
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reactor. 



5. A method of controlling the properties of a product, which comprises computing 
the average operational conditions appropriate for optimum operation of the 

5 . polymer product per lot or batch by using the method of estimating the properties 
of a polymer product according to Claim 1; and then utilizing the same in quality 
control of the polymer product 

6. A method of controlling the properties of a product, which comprises computing 
10 the average operational conditions appropriate for optimum operation of the 

polymer product per lot or batch by using the method of monitoring the state of 
process according to Claim 3; and then utilizing the same in quality control of the 
polymer product 

15 7. The method of controlling the properties of a product according to Claim 5 or 6, 
wherein said process is a process which uses a single polymer reactor or a multi- 
stage polymer reactor. 

8. The method of monitoring the state of process according to Claim 3 or 4, wherein 
20 said process is a general chemical process. 

9. The method of controlling the properties of a product according to Claim 5 or 6, 
wherein said process is a general chemical process. 
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